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Abstract 

We apply a new method based upon thermofield dynamics (TFD) to study entanglement of 
finite-spin systems with non-competitive external fields for both equilibrium and non-equilibrium 
cases. For the equilibrium finite-spin systems, the temperature dependence of the extended density 
matrices is derived using this method, and the effect of non-competitive external field is demon¬ 
strated. For the non-equilibrium finite-spin systems, the time dependence of the extended density 
matrices and the extended entanglement entropies is derived in accordance with von Noumann 
equation, and the dissipative dynamics of the finite-spin systems is argued. Consequently, the 
applicability of the TFD-based method to describe entanglement is confirmed in both equilibrium 
and non-equilibrium cases with the external fieds. 


1 Introduction 

Generally, entanglement states show a quantum-mechanically complicated behavior. A new method 
to analyze quantum entanglement using thermofield dynamics (TFD) jj| [5, 3J has been proposed 
in Ref. j4j, and the capacity has been confirmed in Ref. j5]. In this new treatment of quantum 
entanglement with TFD, an extended density matrix has been defined on the double Hilbert space 
(ordinary and tilde Hilbert spaces), and examined for some simple cases [1]. The new TFD-based 
method allows the entanglement states to be easily understood, because the intrinsic elements caused 
by quantum entanglement can be extracted from the extended density matrix in this formulation. 
Consequently, it has been found that the intrinsic quantum entanglement can be distinguished from 
the thermal fluctuations included in the definition of the ordinary entanglement at finite temperatures. 
Based on the analysis presented in Ref. [4], it has been argued that the general TFD formulation of the 
extended density matrix is applicable not only to equilibrium states but also to non-equilibrium states. 
Furthermore, the extended entanglement entropy has been introduced using the extended density 
matrix, and compared with the traditional entanglement entropy for the case of non-equilibrium 
spin system without external field mm- However, analysis of the entanglement entropies of non¬ 
equilibrium systems with external fields has not been conducted in Refs. sis and, therefore, the 
examination of the entanglement entropies of non-equilibrium systems with external fields using TFD 
is of current interest. In the present communication, therefore, the extended density matrices and the 
entanglement entropies of finite-spin systems with and without the non-competitive external fields are 
exhaustively investigated in both the equilibrium and non-equilibrium cases, based upon the general 
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TFD algorithm. 

The rest of the paper is organized as follows. In the next section, we introduce the extended density 
matrix in the double Hilbert space of the TFD, and review the equilibrium systems with external fields. 
In Sec. 3, we obtain the the extended density matrices and the the extended entanglement entropies 
of non-equilibrium spin systems with and without the non-competitive external fields, and discuss the 
numerical results of the time dependence. Finally, conclusions are described in Sec. 4. 


2 The extended density matrix metod and the example of 
equilibrium finite-spin systems 


First of all, we give a short introduction of the extended density matrix, p, in the double Hilbert space 
of the TFD. p has been defined in Ref. 4] as follows: 

p ■= i^x^i, w ; = p l/ 2 Y I s ’ s ) = p l/ 2 Y l s ^’ fi) 

S S 

using the ordinary density matrix, p, where {| s)} is the orthogonal complete set in the original 
Hilbert space and (|s)} is the same set in the tilde Hilbert space of the TFD [7][9j. If entanglement 
subsystems A and B are being examined, each of |s) and |s) states are represented as the direct 
products |sa,Sb) = |sa)|sb) and |sa,sb) = |sa)|sb}> respectively. Moreover, it has been proven in 
Ref. [4] that if the square-root density matrix p 1//2 is written using the matrix elements a aAt p B:Cl ' At p^ 
as 

p ^ 'y ' a aA,/3B,ct A ,/3 B l a A) ^b}(o:A)/^bI) (2) 

ot A ,pB,a' A ,P' B 

the renormalized extended density matrix pa is then expressed as 


p A := Tr B p := Y (tb,7bIp|7b,7b) 

7b, 7b 

'y 1 I^Aj /3 a) (c^Aj Pa\ 

aA,/3A,a A ,/3 A 

Y K A ,p A ,cc' A ,p’ A (|o:a)(Q:a|) (|^a)<^|) , 
a A >0 A , a ' A iP' A 


where 

ba A ,p A ,a' A ,P' A = a “A,7B,^A,7B a Q!^,7B,/3^,7B' 

7b,7b 

Let us reexamine the S = 1/2 spin system described by the Hamiltonian 


(3) 

(4) 


H:=-JS a -S b -v b H(S z a + S£), (5) 

incorporating the spin operators, S a = ( SS j[) and S b = (fig, S^, S'g), of the subsystems, A and 
B, respectively, where H is the magnitude of external field. The state, |s), of the total system is then de¬ 
noted by the direct product, |s) = |sa, sb) = |sa)|sb)- Using the base, {| + +), | H—}, | —h), | }}, 
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the matrix form of the Hamiltonian (0 is then expressed as 


^ ^ ' ^SAjSB.s^Sg l s Ai s b)( s Ai s 


' J I 
bI 


sa,sb,s a ,s b 


J 


J 


— 4 — VbH j I + +)(+ + I + — + PbH j I-}(-I 

+ — (| + -)(+ - | + I - +}(^ + 1)- 2 (1 + -)(- + I + I ~ +)(+ ~ I) • 


( 6 ) 


As can be seen from Eq. ©, the spin inversion symmetry of this system is broken by external field, 
H , in the Hamiltonian, "H. 

For the equilibrium states in terms of T~L expressed in Eq. © , the ordinary density matrix, 
p eq := e~P n /Z(f3), of this system can be obtained as 

p-k /4 


Pe q — 


m 


e K /2+h.\ + + ^ + + | + e JST/2-fc| _ _ 


+ COSh y (I + -)<+ - | + | _ +)(_ + |) 
+ sinh|(|-+}(+-I + I + -X- + I) 


where j3 is the inverse temperature, the partition function, Z(f3) := Tre /3W , is given by 


Z(/3) = 2e K / 4 ^e K / 2 cosh h + cosh y^j 


(7) 


( 8 ) 


K := (3J and h := j3psH , respectively. For p eq in Eq. 0, the square-root density matrix, p e q 1 ^ 2 , is 
then expressed as 




i Je h ( e h + e~ K + 1) + 1 

1 

\/e h ( e h + e~ K + 1) + 1 

P h / 2 / 1 


+ +X++I 




+ 


2 y y/e 2h + e h + 1 + e h ~ K \Je h 4- e K ( e 2h + e h + 1) 

x(l + -X+-| + l-+X- + l) 


oh /2 


y Ve 2h + e h + 1 + e h K s/e h + e K ( e 2h + e h + 1) 

x(| + -}(- + | + |-+}(+-|). 


(9) 


It is of course that p e q 1 / 2 in Eq. 0 satisfies (peq 1 ^ 2 ) 2 = Peq- According to formulae 0 and 0, we 
then obtain the extended density matrix p ^ := Ti'B/5 eq := (sb, SB|P eq l s B, s' B ) of the spin A, 


p7 = &d q il+)(+ll+)(+l + &d q 2 l-X-ll-X-l 
+ Ki (l+XHI+X-l + I-X+IFX+I) 

+ & qe (l+X+ll-X-l + RHI+X+l) . (10) 
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for the equilibrium system, where the matrix elements and 6® q are respectively obtained 

as analytic functions 


and 


e h (^e h+K +{e K / 2 + l) 2 ) 
dl 4 ( e h + e K ( e 2h + e h + 1)) ’ 

(ID 

(e*/ a + l)V + 4e* 

d2 4 ( e h + e K ( e 2h + e h + 1)) ’ 

(12) 

eq cosh(f) (e^ + 1) e( 2ft + 3A ')/ 4 
bct e h + e K (e 2h + e h + 1) ’ 

(13) 

sinh 2 (f) 

jeq _ V 4 / 

qc 2 cosh ( y) (cosh ( h ) + 1) + 2 sinh (A) cosh ( h ) 

(14) 


of /3 and H , and correspond to the two diagonal elements (dl and d2), the classical and thermal 
fluctuations (cf) and the quantum entanglements (qe) of p respectively. The numerical results of 
these matrix elements at psH/J = 0 and 0.3 have been obtained in Ref. [1], and the breaking of 
the spin inversion symmetry caused by the external field, H > 0, has been argued, i.e. &^ q ^ 6^ a t 

H > 0, and b^j = b e 'y 2 at H = 0. The parameter, 6® q , in Eq. m expresses the quantum entanglement 

effect. This quantum fluctuation is essential to quantum systems, and it has been used as an order 
parameter of 2D quantum systems DUE]- As can be seen from Eq. GOD, only the intrinsic quantum 
entanglement is extracted clearly in the TFD formulation. In particular, it can be understood that 
the entangled state of the system emerges through a single product, such as |+)(+||—)(—|, in pa■ Let 
us now parenthetically touch upon the physical significance of ^di>^d 2 ^cf an d ^qe- As can evidentlly 
be seen from Eqs. gii. m, the high temperature limits of 6^, 6^ q , an( i ^qe are 

lim bj'i = lim bVk = lim 6®? = — and lim 6® q = 0, (15) 

/3-S-0 /3-S-0 /3—s-0 ct 2 /3-s-O q 


respectively, for any value of H. On the other hand, the low temperature limits of &)) q , f\y 2 . 6® q and &® q 
are 


lim foj'} = 1, lim b e 3 = lim b e 3 = 0 and 

/3->oo di /3-s-oo /3—»oo cl 


lim 6® q = 0, 
/3—too 


respectively, for H >0, and 


lim 

0—>oo 


b eq 

a dl 


lim 

/3—>-oo 


b eq ~ 

u d2 ~ 
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12 ’ 


lim 6®? 
j3 —>-oo 


— and lim 6® q = —, 
3 /3^oo q ® 12 


(16) 


(17) 


respectively, for H = 0. These results are summarized in Fig. G1 and clarify that a transition 
from broken phase (H > 0) to unbroken phase (H = 0) of spin inversion symmetry happen in the 
equilibrium entanglement finite-spin system. 
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Figure 1: Temperature dependence of the elements b^, b^ 2 . b^j. and in Eqs. (fill) . ~, (flTl) . Parts 
(a), (b), (c), and (d) show cases with ^bH/J = 0.1, 0.01,0.001 and 0, respectively. The dashed lines 
in all the parts represent the asymptotes of the 6^, b^, and curves, and the asymptote of the b ®a 
curve is the holizontal axis. In parts (a), ~, (c), at /3 —>■ oo, and b ®^ converge to the same 

value 0, on the other hand, converges to 1. Part (d) shows that = 6^, b ^, and &®a converge to 
Y : ' 2 , y and jh, respectively, at /3 —> 00 . 
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3 Non-equilibrium finite-spin systems with non-competitive 
external fields 

Let us consider a non-equilibrium system with dissipation, which is described by the Hamiltonian of 
Eq. (O. We assume that the time dependence of the ordinary density matrix, p(f), of this system is 
given by the dissipative von Neumann equation EJlz], which is described as 


iH §t P ^ = ~ 6 ~ Peq ^ ’ 

with e being a dissipation parameter. The solution of Eq. m is then expressed as 

p(t) = e~ et U ] {t)p 0 U{t) + (1 - e~ et )p eq , 

for any initial condition, p(0) = po , where the unitary operator, U(t) := e znt t h , denotes 
U(t) = eW (exp | + +)(+ + | 

+ exp(^^®) l -- ) (--| 

+ COS y (I + -}<+ - | + | - +}<- + I) 

-isiny(| -+)(+- | + I + -X- + I) 


(18) 


(19) 


( 20 ) 


and ui := J/H. The explicit expression of p(t) in Eq. (1191) is complicated for any initial condition, so 
hereafter, we confine ourselves to the initial condition, pg = | H—}(-|— |. The insertion of Eqs. Q) 
and ®, along with the initial condition, into Eq. OH), brings forth 


P(t) = 


+ 


+ 


2 (e et - 1) e K+2h 
e K ( e 2h + e h + 1) + e h 
2 (e et - 1) e K 
e K ( e 2h + e h + 1) + e h 

(e et — 1) cosh (y) 
cosh (y) + cosh (h)e K /‘ 2 

{e et — 1) cosh (y) 


+ 


cosh (-|-) + cosh (h)e K / 2 

e h ( e et — 1) {e K — l) 
e K (e 2h - + e h + 1) + e h 

' e h ( e et - 1) ( e K - l) 
e K ( e 2h + e h + 1) + e 


+ +)(++ I 
I-)(-I 

+ coswf + 1^ | H-)(d— 

— cos cat + 1 ] | —1-)(—h 


— isinwtj | -|—)(—b 
h + isinu>t ) | - +)(+ - 


( 21 ) 


where h := pb/3H/(u}H). 

The ordinary entanglement (von Neumann) entropy, S := —fceTrA [pAlogpA] and pa ■= Tib p(t), 
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are then paraphrased into 


* = -T e 


x log 


fc B ( ( ((e h +2) e K + e h ) e u + ( e 2h — l) e 


K 


(e h +1 


,2h 


1 )' 


— COS Lot 


(( e h + 2) e K + e h ) e t£ + (e 2/l — l) e 
( e h + e 2h + 1) e K + e h 


— cosuit 


(2e h+K + e K + l) e h+t£ — ( e 2h — l) e 


K 


cos cot 


x log 


(e h + e 2h + 1) e K + e h 

/ (2e h+K + e K + 1) e h+te - ( e 2h - l) e 
l (e h + e 2h + 1) e K H 


K 


COS LOt 


( 22 ) 


It is also worth mentioning that S in Eq. (l22l) is proportional to an entanglement, E(C), which is a 
function of the “concurrence” [H] , 

P -tt 

C = 


\/2 (e K ( e 2h + e h + 1) + e h ) 
x ((- (2^+* + + 1) e et+h ~{{e h + e 2h + l) e K + e h ) cos (too) + ( e 2h - l) e K ) 

x (_ e e* (( e * + 2 ) e x + e h ) + ((e h + e 2h + l) + e h ) cos(tw) - (e 2fl - l) e*)) :1/2 . 


(23) 


The time dependence of S and C is displayed in Fig. [3] (in units of fee = 1)- In the dissipative system, 
at t —> oo, S and C converge to the values, 


2e K + e K+h + e h 
2 (e K ( e 2h + e h + 1) + e h ) 
e h ( e K + 2e K+h + l) 


log 


log 


2e K + e K+h + e h 
2{e K (e 2h + e h + 1) + e h ) 
e h (e K + 2e K+h + l) ^ 


and 


2 (e K (e 2h + e h + 1) + e h ) l 2 (e K ( e 2h + e h + 1) + e h ) 


\j2e K+2h (2 sinh/\ + (2 cosh (h) + 3) cosh A' + 2 cosh(/i)(sinh K + 1) + 1) 


(24) 


(25) 


e K (e 2h + e h + 1) + e h 

respectively, so it is possible to argue that S and C include not only the contribution of the quantum 
entanglement, but also the contribution of the classical and thermal fluctuations in the equilibrium 
case. However, this fact is not clarified in Eqs. (12211 and (1231 . 

We are now ready to study the extended density matrix in the TFD double Hilbert space. In terms 
of the expression similar to Eq. m, we are then led to the renormalized extended density matrix, 

~neq 

Pa > as 


PA eq = 6STI+K+II+K+I + ^ q l-)(-ll-)(-l 

+ KP (l+X-ll+X-l + l-X+ll-X+l) 

+ foqS q (l+x+ll-x-l + l-x-ll+x+l), (26) 

where the matrix elements, ^cf q anc ^ ^qe q ; are respectively expressed as analytic functions of 

H , e, t and /3. However, their expressions are so complicated that we show their numerical results of 
several cases in Fig. [2j 
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Figure 2: Time dependence of matrix elements, 5)[® q , b^ 1 , q and b q l q , in dissipative and non- 
dissipative systems with scaled temperature, (/3J) _1 = 0.7. Parts (a), (b) and (c) show cases with a 
scaled dissipation rate of e/w = 1, 0.1 and 0.01, respectively. The dotdashed, dotted, dotdotdashed 
and dashed lines in parts (a), (b) and (c) represent the asymptotes of the ^di q , q ,and fe"| q 
curves, respectively. In part (d) which is the non-dissipation case (e = 0), 5“ q vanishes identically. 


The extended entanglement entropy, S := — /c b Ti-a [/5“ q log p” q ], in the present case is then re¬ 
duced to 

S = S d + S qe , (27) 

where 



and 


S qe := -2fc B 6qI q log6“ q , (29) 

respectively. In Eqs. 03, (EH) and m, the expressions of S, the classical and thermal fluctuation 
parts, Sd, and the quantum entanglement part, S qe , also comprise analytic functions of t, /3, e and to, 
respectively, however, their expressions are quite complicated. So, we show the numerical behaviour 
of C, S, S, S qe and 6“| q for a few cases in Figs. [31(a) — (c) (in units of = 1). These figures 
exhibit that, at t — > oo, S, S qe and foq| q converge to the values, 0.277 • • • , 0.150 • • ■ and 0.019 ■ • •, 
respectively. The traditional entanglement entropy, S, consequently becomes larger than the extended 
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Figure 3: Time dependence of C, S, S and S qe , along with parameter &q| q , in dissipative and non- 
dissipative systems with (/3J) _1 = 0.7 and ^bH/J = 0.3. Parts (a), (b) and (c) show cases with a 
scaled dissipation rate of e/w = 1, 0.1 and 0.01, respectively. The dotdashed, dotted, and dashed 
lines in parts (a), (b) and (c) represent the asymptotes of the C, S and S qe curves, respectively. In 
(a), ~, (c), 6“ q converge to the value 0.019 • • •. In part (d) which is the non-dissipation case (e = 0), 
all the curves have the same phase. 


entanglement entropies, S and S qe , at t —> oo. These results suggest that these numerical behaviour 
of C, S, S, S qe and 6“ q is resemblant to the one in the case of no external field [5]. 

For non-dissipative systems, S in Eq. (1271) and i§ qe in Eq. (1^1) reduce to 


S 


-/SB ( cos 4 y • log (< 


cut 




cos 


. 4 ut , 

sin — • log 
2 6 


1 


+ - sin Lot ■ log 


sin cot 



(30) 


and 


Sqe = -y sin2 Ut • log , (31) 

respectively, at e = 0. As can transparently be visualized from Eqs. (l30l) and (1311) . the dependence of 
S and S qe on both J3 and H disappear at e = 0 and this result assorts with the case of no external 
field. The time dependence of S and S qe at e = 0 is exhibited in Fig. [2|d) (in units of fee = 1). It 
is apparent in this figure that all the curves (C, S, S , 5)^ and b qe ) showing the entanglement have 
the same phase, however, their amplitudes differ. This result is also similar to the case of no external 

field [5]- 
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4 Conclusions 


In this communication, we have examined the extended density matrices and the extended entangle¬ 
ment entropies of finite-spin systems with and without non-competitive external fields in both the 
equilibrium and non-equilibrium cases, based upon the TFD formulation. For the equilibrium system, 
the extended density matrix elements are derived using TFD, and the low and high temperature lim¬ 
its are calculated explicitly. The results show that the effect of the non-competitive external field is 
significant for the temperature dependence of the extended density matrix. For the non-equilibrium 
systems, based upon the extended density matrix, C, S, S, S qe and &q| q , are computed using TFD, 
and summarized in Figs. [2] and [3j From these results, we can deduce the following: 

1) The each shapes of the curves of b r ( y 2 q and b q l q in Fig. [2]^b) are similar to those of Fig. 4 in 

Ref. ^ , except that the curve of q is compartively wavy. 

2) The each behaviours of &^® q , >^cf q an( i ^qe q ai Figs. (21a), (c) and (d) are unparalleled. 

3) The each shapes of the curves of C, S, S, S qe and 6q® q in Figs. (31a), (b), (c) and (d) are also 

similar to those of Figs. 1 (a), (b), (c) and (d) in Ref. [5]. 

As a consequence, the non-competitive external field has been found to have little effect on C, S , S, S qe 
and 6q® q , and modest effect on So another thing to keep in mind is to examine another spin 

system with competitive external field using the TFD formulation. However, the dissipative dynamics 
of entanglement in frustrated spin systems is more complicated, and is under consideration. 

We have anyhow succeeded in elucidating the extended density matrices and the entanglement 
entropies of finite-spin systems with and without the non-competitive external fields which are con¬ 
structed in both the equilibrium and non-equilibrium cases, based upon the general TFD algorithm. 
Let us conclude by emphasizing that the new TFD-based method enables us to clearly distinguish 
between the various states of quantum systems. 
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